Raster

Layer functionalities

Filters
Description

Filtering is a process by which we can enhance images. gvSIG can filter images through a variety of filtering methods. In the upper left part of the Filter dialog, the filters are grouped by type (1). By double-clicking one of the filters or by clicking on the "Add Filter" button on the bottom left, the filter will be added to the list of filters in the lower left part of the Filter dialog. All filters in the filter list will be applied in the preview. If you want to remove a filter from the list, you can either double-click on the filter or click on the "Delete filter" button. The filters in the list will be applied to the image in the order that they appear. Keep in mind that the order in which the filters are applied will affect the result, and changing the order of the filters may change the output.

In the middle of the dialog window are the controls of the selected filter (2). When changing the controls of one of the filters from the filter list, the results will be directly shown in the preview window. Below the middle part of the dialog you can change the name of the output layer that will be generated when clicking "Apply" or "Close".

On the right side of the dialog you can preview the outcome of the filters (3). (See documentation on "Preview tool"). In the lower right part you can select whether you want to display the filters over the selected layer or save the filtered image as a new layer (4).

The button "Apply" will apply the changes according to the entered parameters, keeping the Filter dialog open. The "Close" button will apply the changes and close the Filter dialog. The "Cancel" button will close the Filter dialog without applying any filters.

All filters in the filter list can be activated or de-activated through the "Active" checkbox. This checkbox is usually located in the upper part of the filter control panel.

filtrado.img/es/Filtros.png

Configuration panel for the image filters

Generate a new layer or apply to current layer

The number of applied filters will affect the time that it will take to draw the layer. If you choose to apply the filters to the current layer, the drawing and re-drawing of the layer may slow down while the filters are applied. If the filter results are saved as a new layer, the filtering process has to be done only once so that the next time the layer is drawn, it will not be slowed down by the filtering. Therefore, it is generally recommended to save the output to a new layer if possible. There are cases though in which it is not recommended to generate a new layer. For example, if you have a large orthophoto and you only want to change the brightness a little, it could take more time to save the output as a new layer. If the brightness filter is applied over the current view, the area on which the filter is applied is much smaller which makes the drawing faster. It is up to the user to decide whether it is better to create a new layer or display the filters on the view of the current layer.


Enhancements

The brightness filter changes the brightness value of the layer. You can increase or decrease the brightness by moving the position of the sliding bar or by entering the value directly in the text box and press enter.

filtrado.img/es/Brillo.png

Brightness filter

The contrast filter changes the contrast value of the layer. You can increase or decrease the contrast by moving the position of the sliding bar or by entering the value directly in the text box and press enter.

filtrado.img/es/Contraste.png

Contrast filter


Spatial filters

With this type of filter, graphical transformations like smoothing, edge detection, sharpening etc. are applied to the image.

The following filter types can be applied:

MEDIAN FILTER

The median filter applies a kernel of a certain size, which is determined by the user through the sliding bar labeled Window side.

The median filter is normally used to smoothen and to reduce noise in an image, by moving a kernel of N x N number of pixels over the image and evaluating each central pixel, replacing its value with the median of its neighboring pixels. Compared to the Mean filter, the advantage of the Median filter is that the final pixel value is a value that actually occurs in the image and not an average.

filtrado.img/es/median-es.png

Median filter

MEAN FILTER

The mean filter applies a kernel of a certain size, which is determined by the user through the sliding bar labeled Window side.

The filter replaces the value of the central pixel with the mean value of the surrounding pixels. Each value of the kernel would be one and the divider would be the total number of elements in the kernel (i.e. a kernel of 3 x 3 would replace the value of the central pixel by the average value of the nine pixels covered by the kernel).

filtrado.img/es/Media.png

Mean filter

LOW PASS FILTER (smoothing filter)

The low pass filter applies a kernel of a certain size, which is determined by the user through the sliding bar labeled Window side.

Using a low pass filter tends to retain the low frequency information within an image while reducing the high frequency information.

filtrado.img/es/PasoBajo.png

Low pass filter

SHARPENING FILTER

By moving the slider to change the sharpness (values from 1-100), the contrast of an image can be changed. The results can be evaluated in the preview window. With a higher contrast, details in the image can be accentuated but the noise will also increase.

filtrado.img/es/Sharpen.png

Sharpening filter

GAUSS FILTER

The Gauss filter applies a kernel of a certain size, which is determined by the user through the sliding bar labeled Window side.

The maximum value appears in the central pixel and gradually decreases for pixels that are further away from the central pixel.

filtrado.img/es/Gauss.png

Gauss filter

CUSTOM FILTER

This is a kernel of 5 x 5 or 3 x 3, for which the values can be introduced by the user. After multiplying the pixel values with the kernel values, the result will be divided by the number specified in the Divisor textbox.

filtrado.img/es/Personalizado.png

Custom filter

MODE FILTER

The mode filter applies a kernel of a certain size, which is determined by the user through the sliding bar labeled Window side.

This filter takes the value that occurs most in the surrounding pixels and assigns it to the central pixel.

filtrado.img/es/Moda_es.png

Moda filter


Colour adjustment

Adjustment of RGB values

It is possible to change the balance between Red, Green and Blue in an image if needed. To do this, move the sliding bar to increase or decrease the values or enter the value directly in the text box next to the sliding bar. Ticking the "Brightness" check box ensures that the brightness level of the pixels will be maintained while the RGB values are changed.

filtrado.img/es/RGB.png

RGB balance filter

Adjustment of CMY values

It is possible to change the balance of Cyan, Magenta and Yellow in an image if needed. To do this, move the sliding bar to increase or decrease the values or enter the value directly in the text box next to the sliding bar. Ticking the "Brightness" check box ensures that the brightness level of the pixels will be maintained while the CMY values are changed.

filtrado.img/es/CMY.png

CMY balance filter

Adjustment of HBS values

It is possible to change the balance of Hue, Brightness and Saturation in an image if needed. To do this, move the sliding bar to increase or decrease the values or enter the value directly in the text box next to the sliding bar.

filtrado.img/es/HSL.png

HBS balance filter


Edge detection

These filters attempt, through the use of kernels, to detect edges in the image and change the image so that these edges are enhanced, while the rest of the image is grayed out.

filtrado.img/es/filtro_bordes_es.png

Filter dialog. Edge detection

There are four edge detection filters, all with the same interface and options, in which the user chooses a threshold in the range 0-255, and the possibility compare the results by ticking the compare check box:

filtrado.img/es/Bordes.png

Sobel filter example

SOBEL

The Sobel filter detects the horizontal and vertical edges separately on a grayscale image. Colour images are converted to RGB gradations. The result is a transparent image with black lines and some remains of colour.

ROBERTS

The Roberts filter is suitable for detecting diagonal edges. It offers good performance in terms of location. The major drawback of this filter is its extreme sensitivity to noise and therefore has poor detection qualities.

PREWITT

The Prewit filter detects edges in all directions as it consists of 8 kernels that are applied over the image pixel by pixel.

FREI-CHEN

The Frei-Chen filter processes the neighbouring pixels as a function of their distance from the pixel that is being evaluated. The result is that edges in all directions are detected.


Masks

Transparent area

With this functionality it is possible to set the transparency level of a Region of Interest (ROI). The region of interest must have been defined previously. If the layer does not have a region of interest, the following message will appear: "A Region of Interest (ROI) must be defined for this layer to apply this filter. Please go to the dialog Area of Interest and select at least one ROI." If there are already one or more ROI associated with the layer, the message will not appear. Instead, a list of ROI will be shown, from which you can select one or more by ticking the corresponding check box. Then, adjust the level of transparency with the slide bar or by entering the value directly in the text box next to the slider. Ticking the check box labeled as "Inverse" will result in the opposite effect; all of the image except for the ROI will be set to the specified transparency level.

filtrado.img/es/RegionesTransparentes.png

Transparent area filter

Mask

With this functionality it is possible to cut out a Region of Interest (ROI) that has been previously defined for the layer by assigning a fixed user-specified value to the rest of the image outside the ROI. If the layer does not have a region of interest, the following message will appear: "A Region of Interest (ROI) must be defined for this layer to apply this filter. Please go to the dialog Area of Interest and select at least one ROI." If there are already one or more ROI associated with the layer, the message will not appear. Instead, a list of ROI will be shown, from which you can select one or more by ticking the corresponding check box. Then, select the value to be assigned to the pixels outside the ROI by typing a number in the "value" text box. The default value is -99,999. Ticking the check box labeled as "Inverse" will result in the opposite effect; the ROI will be assigned the specified value while the rest of the image values are maintained.

filtrado.img/es/Mascara.png

Mask filter


Histogram
Description

To launch the histogram dialog window, use the drop-down toolbar selecting the "Raster Layer" button on the left and "Histogram" in the drop-down button on the right. Make sure that the text box that displays the current layer is set to the name of the raster layer for which you want to see the histogram.

histograma.img/es/barraHistograma.png

Histogram icon

The Histogram dialog shows a histogram of the statistical distribution of pixel values in the current view. This information is often useful when you are trying to color balance an image. In the middle of the dialog you will see the graph on which you can right-click to show a context menu with general options for this kind of graphics.

histograma.img/es/histograma.png

Histogram dialog window

In the upper part of the dialog (1) are the controls to configure the histogram:

1. Type of histogram

There are three types: "Normal", "Accumulated" and "Logarithmic".

2. Data source

With this option you can select the data source for the histogram:

Current view (R,G,B):

With this option, the pixel values that are displayed in the current view of gvSIG will be used for the histogram. Therefore, the band selector shows only the R, G, and B values which are the visual bands. Every band will appear in its corresponding colour in the graph (red for R, green for G and blue for B). This is the default option when the histogram dialog is opened.

Complete histogram:

With this option, the histogram for the whole raster layer is calculated. Because of the amount of time that it would take to calculate the histogram for large images, the histogram is only calculated once and saved with a .rmf extension in the directory in which the image is stored. After the first time, the histogram for the same layer can be displayed much faster. (Keep in mind that if you delete the .rmf file that is stored with the image, you will lose its histogram information.)

3. Band selection

Apart from identifying to which band each histogram corresponds through its colour (in case of the current view Data Source) you can also identify the band by hovering the mouse over a point in the graph. The tooltip displays the band name and the value of the point.

4. Context menu and manipulation of the graph (2)

Zoom operations

We can zoom in and out of the graph using the mouse.

You can also zoom in and out using the context menu.

Context menu

When you right-click on any part of the graph, the context menu is shown with the following options:

histograma.img/es/menuContextual.png

Histogram options. Context menu

histograma.img/es/propiedadesGrafico.png

Histogram properties

5. Statistics (4)

The controls that appear under the graph allow the user to restrict the range of values (X axis of the histogram) on which the histogram is based. The default setting is the complete range so that, for example in a Byte data type image, the statistics are calculated for all the pixel values from 0 to 255. You can enter the values directly in the text boxes or use the + and – controls next to the text boxes. You can also slide the triangles over the sliding bar to select the range of values.

histograma.img/es/barraEstadisticas.png

Sliding bar with pixel ranges

In this table, the statistics that correspond to the selected range of pixel values are shown in the text boxes. Each row of the table corresponds to one raster band as displayed in the histogram. The columns that are shown are:

6. Export the table (3)

You can export the table through the option "Save as DBF". The data contained in this table are the values of the current histogram. After creating the DBF table, it can be used as any other table in gvSIG.

histograma.img/es/tabla.png

Resulting DBF table


Preferences

The Raster section of the Preferences dialog contains the option "Number of classes" where you can set the number of intervals in which the histogram is divided when the data type of the image is not Byte. For Byte images, this value is 256. In the preferences dialog, the default value of this option is 64 but you can choose any of the options (32, 64, 128, and 256). The intervals are the parts in which the range of values is divided. For example, if we have a DTM with values between 0 and 1 and there are 64 intervals, each interval will have a range of 1/64.

The number of classes does not only refer to histograms but also to other functionalities that require a division in intervals of value ranges.

histograma.img/es/preferencias.png

Raster Preferences


Radiometric enhancement

The maps that are obtained through digital processing of satellite imagery are useful not only for thematic mapping, but also as a backdrop on which map features can be overlaid. If the visible bands are displayed in a colour composition through the colouring of each band with the corresponding colour gun, it is important that the bands are sufficiently enhanced so that the colours appear more natural. The final display colour depends not only on the direct result of the chosen colour composition but also on the radiometric post-processing. The satellite image map will be more useful as backdrop if the bands are enhanced and displayed in colours that match the natural colours as the human eye perceives them. gvSIG provides the enhancement tools to adjust the colours for each band.

In the following sections the different parts of the dialog are described:

Histograms

The central part shows two graphs (1). The graph on the left is the histogram of the input image. The graph on the right shows the histogram of the output image. The graphs that are presented with a yellow line can be modified with the mouse. When you change the input histogram, the output histogram will be changed accordingly and you can preview the result.

In the upper corners of the input histogram are the maximum and minimum values of the raster displayed. In the lower corners, the maximum and minimum values that are being included in the enhancement are displayed. The percentage of values that are being left out of the histogram appears in parentheses. These values can be modified by grabbing and dragging the dotted vertical lines on the side of the graph. Dragging the left line will modify the minimum value, while dragging the right line will modify the maximum value. (This way, by leaving out the values that are not used in the input image, you can stretch the output values over the whole range of available values, so that the visual quality is improved.)

realces-radiometricos.img/es/levelslice.png

Radiometric Enhancement dialog

Controls

In the lower part of the dialog (2) you will find some controls with the following options:

Type of function:

The enhancements will replace each input value with an output value. This process is done by creating a look-up table which provides the correspondence between a range of input data and a range of output data. To apply this correspondence, a fuction is used. The used function and its parameters are chosen by the user.

Linear enhancement

realces-radiometricos.img/es/linealRaster.png

Linear radiometric enhancement

realces-radiometricos.img/es/density-slicingRaster.png

Piecewise linear enhancement

Non-linear enhancement

The non-linear enhancements have the same approach as the linear enhancements in the sense that each input value is replaced by an output value. The difference lays in the function that is assigned to produce the output values, which is non-linear. The available non-linear functions are logarithmic, exponential and square root. With each function you can modify the curve to smooth or accentuate the enhancement result.

realces-radiometricos.img/es/exponencialRaster.png

Exponential radiometric enhancement

Band

With this option you can specify the raster band to which the enhancements are applied. For a correct balance of the image, it is recommended to enhance each band separately.

Drawing type

With the option drawing type, different types of histograms can be chosen. Filled will draw a filled histogram while Line will only show the contours of the histogram. The colour of the line or fill pattern depends on the selected band. The bands Red, Green, Blue and Gray are displayed in red, green, blue and gray respectively.

Type of histogram

RGB Check box

When check box labelled as RGB is ticked, it is assumed that the image is displayed as RGB with Byte data type and values between 0 and 255. If the checkbox is not ticked, it is assumed that the range of values are Byte data type values between -127 and 128, which will produce significant differences in the display and in the minimum and maximum values that are shown in the bottom of the input graph.

Display enhancement results

In the lower right part of the dialog (3), you can indicate how you want to see the enhancement results; in the current view or saved as a new layer.

Preview

The preview window (4) shows the real-time results of each enhancement that is applied to the image.


Save View as image

The tool for exporting the view as an image can be accessed from the drop-down toolbar by selecting "Export to raster" on the left button and "Save view to georeferenced raster" on the right button. Make sure that the name of the raster layer that you want to export is set as the current layer in the text box.

salvar-a-raster.img/es/menu.PNG

Export to raster. Save view to georeferenced raster

A message will appear to inform that you can use the selection tool to set the area in the view to export.

salvar-a-raster.img/es/confirmacion.PNG

You can begin to select a area on the view

Now, you can select two points in the view to define the rectangle of the area to be exported, by clicking the first point and dragging the mouse towards the second point, then release.

salvar-a-raster.img/es/seleccion.PNG

Selection of a rectangle to define the output image

Then, the Save view to georeferenced raster dialog will appear. If the selected area is too small, the dialog will not appear and a bigger rectangle must be selected.

salvar-a-raster.img/es/captura1.PNG

Save view to georeferenced raster dialog

The upper part of the Save view to georeferenced raster dialog shows the coordinates of the two points that define the selected area in the view. You can edit the coordinates to change the selected area.

In the option box in the central part of the dialog you can choose from three selection methods:

NOTE: To save time and memory the maximum size of output images is limited to 20000 x 20000 pixels. If the intended output image is larger and you click on "Apply", gvSIG will display a message that the parameters must be changed before trying again.

Clicking the "Selection" button will open a file browser dialog where you can specify the output file. Depending on the type of file, the corresponding driver will be loaded (you will notice that the button on the right of the "Selection" button will change). For example, an output file .jp2 will open the properties dialog for Jpeg2000. The formats in which you can save are .TIF, .IMG, .BMP, .PGM, .PPM, .MPL, .RST, .JP2, .JPG, and .PNG. Furthermore but only on Linux kernel 2.4 you can also select ECW.

salvar-a-raster.img/es/captura2.PNG

File browser dialog to save the output image

When you select the output file, the Properties button will be enabled.

For example, for geoTiff the dialog will look like this:

salvar-a-raster.img/es/captura3.PNG

Properties geoTiff

When the output image is selected and the properties set, you can click on "Apply". A progress bar will appear. Depending on the size of the output file, this process may take while. Processing times may vary between a few seconds or several days, so it is important to check the size of the output image in pixels before clicking "Apply". When finished, a screen with statistics will appear that indicates the path of the output image, the disk size, the duration of the process and whether it was compressed. To check the georeferencing of the output image, you can add it to the view as a new layer with transparency.


Clipping layers

(The clipping tool can be accessed from the raster toolbar by selecting "export to raster" from the left drop-down button and "Clipping" from the drop-down button on the right. Make sure that the raster layer that you want to clip is set as the current layer in the text box.)

With the clipping tool, you can create new layers from an existing one. The options are:

Selection of the clipping area

In the "Coordinates" tab of the clipping dialog, there are text boxes to enter coordinates. In the upper part are the values in pixel coordinates and in the lower part the real coordinates. For each item, the two upper text boxes correspond to the coordinates of the upper left corner, while the lower two text boxes correspond to the lower right corner. When changing the pixel coordinates, the real coordinates are re-calculated automatically and vice versa.

There are 3 selection methods that will fill the coordinates automatically. These methods can be activated by clicking the buttons on the bottom of the clipping dialog. From right to left, the buttons are:

recorte-de-capas.img/es/Coordenadas.png

Clipping dialog. Coordinates tab

Modifying the resolution

In the "Spatial resolution" tab of the Clipping dialog, you can modify the resolution of an output image through various interpolation methods. There are two option boxes located on the upper part of this tab:

When modifying the resolution it is necessary to resample and re-assign the pixel values for the output image through an interpolation method. There are four interpolation methods available: Nearest neighbour, Bilinear, inverse distance and B-Spline. The nearest neighbor is the fastest interpolation method, but the results in pixilation of the image and a lower visual quality. The other interpolation methods produce a smoother result.

The button labeled "Restore" returns the initial values of the input image.

recorte-de-capas.img/es/Resolucion.png

Clipping dialog. Spatial resolution tab

Band selection

The "Bands selection" tab of the Clipping dialog displays a table that lists the bands of the input image. When processed, the output image will have the bands in the order as shown in this list. By default, the output image will have the same order of bands as the input image. The order of the bands can be modified through the "Up" and "Down" buttons. The selected row will go up or down one position in the list. The bands can also be omitted from the resulting image by un-checking the corresponding row.

recorte-de-capas.img/es/Bandas.png

Clipping dialog. Bands selection tab

Selecting Options

The "Options" tab of the Clipping dialog presents various options that can be set by the user:

recorte-de-capas.img/es/Opciones.png

Clipping dialog. Options tab


Zoom to raster resolution

You can zoom to raster resolution by right-clicking on the layer in the TOC. In the context menu that appears, click "Zoom to raster resolution".

This will activate a crosshair cursor in the gvSIG view which allows users to perform an action by clicking somewhere in the view. The action in this case is that with every mouse-click, the view will be centered on the point where you clicked. In addition, the view will zoom so that one screen pixel is the same size as a pixel in the current raster layer.


Automatic vectorization

The "automatic vectorization" function can be launched from the raster toolbar by selecting "Raster process" on the left drop-down button and "Vectorization" on the drop-down button on the right. Make sure that the name of the raster layer that you want to vectorize is displayed as current layer in the text box.

vectorizacion-automatica.img/es/barra_vect.png

Vectorization icon

With automatic vectorization, you can generate a vector layer from a raster image using preprocessing to highlight the features of interest.

When launching the Vectorization dialog, the first step is to select the area of the image that you want to vectorize. Keep in mind that the vectorization process may take a long time, so it is recommended to minimize the area (number of pixels) for vectorization. The selection of the area for vectorization can be done in several ways. You can type the coordinates directly; either in pixel coordinates or in the map coordinates. The area can also be selected from the view by clicking the button "Select from the view", after which you can draw an approximate rectangle to define the area. Another selection option is by Region Of Interest (ROI). You can define a ROI here or use a previously defined ROI to set the area for vectorization. In the section "ROI selection" appears a list of available ROI and a checkbox next to each of these to select one or more ROI that you want to use. There are two options to vectorize the ROI: to vectorize the entire area inside the rectangle (bounding box) that covers all the selected ROI, or vectorize only the areas inside the ROI while considering the values outside the ROI as NoData values, excluding them from the calculations.

Finally you can select the scale of the image to preprocess. This is useful because a higher resolution of the preprocessed image will result in a higher precision for the resulting vector layer. You can define this with the drop-down text box labelled "Output Scale". By default, the resolution will be the same as the input image.

When moving on to the next step of the wizard, the process of cutting the image for preprocessing is started. A progress bar appears with the warning that this operation could take a few minutes. The resulting image cut is saved in the temporary folder of gvSIG.

vectorizacion-automatica.img/es/vect_solapa_seleccion.png

Vectorization. Selection of the area to vectorize

There are two methods to preprocess a raster image to vectorize. The first is by creating a limited number of grayscale levels from the original image. The image will be converted to grayscale using one single band or a combination of bands (use the drop-down button labelled "Bands"). For the conversion to grayscale, a posterization process is used to reduce the number of different values. (By default, the image is reduced into 2 levels only: black and white.) For this process you can control the threshold on which the values are passing from black to white and vice versa. This can be done by moving the "Treshold" slider while you can see a preview of the result. (The Treshold slider is only available when there are 2 levels; when there are intermediate grayscale levels, the slider is disabled.) In addition to the posterization threshold, you can apply a mode filter or a noise filter to smoothen the result.

vectorizacion-automatica.img/es/vect_solapa_grises.png

Vectorization. Conversion to grayscale

The second preprocessing method is useful to vectorize contour lines and can be applied to data types other than byte. With this method you can define intervals between each contour line to be vectorized. You can specify the number of intervals in which you want to divide the raster, or indicate the size of each interval. The cuts that have been selected will be shown on a graph that represents the histogram of the image. On this graph, you can modify the distance between cuts, or add or remove some of them using the mouse. It is also possible to modify the distance between cuts in numeric format using the table on the right of the histogram. Each entry in the table represents a cut with the corresponding value. This type of preprocessing is used for digital elevation models (for example .adf or .asc images).

When moving on to the last step of the vectorization wizard, the preprocessed image is generated with the specified values, and saved in the temporary directory of gvSIG.

vectorizacion-automatica.img/es/vect_solapa_tramos.png

Vectorization. Define intervals for vectorization in case of Digital Elevation Models

The last step is to select the method for generation of vectors. There are two methods: contour and potrace, that can be selected from the drop-down button after which a panel appears with settings that are specific for the method. The first method is the simplest and does not have any options. This method will trace the vectors in straight sections going through the pixel centers. This generates a network of vectors based on very small straight sections. The potrace method uses the potrace library for vectorization. The available options for this method are those that the potrace library provides and they are used to define the precision of the tracing of the curves: number of points for each curve, threshold, optimization, etc.

vectorizacion-automatica.img/es/vect_solapa_options.png

Vectorization. Options for vector generation

When clicking on "Apply" or "Accept", the process of vectorization will start after which you will be prompted whether to display the generated layer in the TOC.


Analysis view

The "Analysis view" can be launched from the raster toolbar by selecting "Raster layer" from the left drop-down button and "Analysis View" on the drop-down button on the right. Make sure that the name of the raster layer that you want to analyze is displayed as current layer in the text box.

vista-de-analisis.img/es/barraVistaAnalisis.png

Analysis View icon

With this functionality you can zoom in on the current raster layer with 3 different zoom levels:

  1. On the left of the view, the layer is added to the locator map of gvSIG. This provides a general view of the layer, and you can zoom into the locator map by clicking and dragging, thus drawing a red rectangle. The area inside the red rectangle will be displayed in the view.
  2. The view itself is the second zoom level which functions independently, and the zoom variations that are performed on the view will be applied to the locator map as well so that it keeps being centered on the correct area.
  3. When launching the Analysis view, a small floating window labelled "Cursor Zoom" appears in the upper right corner of gvSIG. This window has the highest zoom level. The zoom level is fixed and always centered on the mouse point. By moving the mouse over the gvSIG view, you will see the contents change.

You can change the relation between the zoom level of this floating window and the gvSIG view. This is done by right-clicking on the floating window and selecting one of the values that are shown in the drop-down menu that appears. The available options are x4, x8, x16, and x32. This means that the pixels in the floating window will be 4, 8, 16, or 32 times bigger than the original.

The floating window also shows the RGB values of the pixel on which the cursor is currently located. The text colour of the RGB values as well as the colour of the central cross (red by default) can be changed by right-clicking on the floating window and choosing the option from the drop-down menu.

Keep in mind that, to see the effects in the floating window while moving the mouse over the view, the view must be active. If it is not active, just click on the view. When the cursor is outside the view, the content of the floating window appears black.

vista-de-analisis.img/es/VistaAnalisis.png

Screenshot with the different elements of Analysis View

There can only be one Analysis view open at any time in gvSIG. Therefore, the button "Analysis View" is re-labelled as "Close Analysis View" when the Analysis view is already open, so that it can be closed before re-opening.


Geographic Transformations

Geolocation

The Geolocation tool provides the ability to change the related transformation that is applied to a raster in its display. A raster could have coordinates that place it in a geographical position. This geographic location is only a change of position, scale and rotation of the image relative to an original position. These changes in position, scale and rotation are those that can be changed with this tool. This can be done numerically and by visual approximation.

This tool does not provide capabilities to deal with complex projections and georeferencing involving the need for resampling.

To launch the geolocation dialogue of the layer, the drop-down toolbar is used by selecting "Geographical Transformations" on the left button and "Geolocation" on the dropdown button on the right. Make sure that the text shows the name of the layer to which we like to assign the transformation.

geolocalizacion.img/es/barraGeolocalizacion.png

Geolocation Tool

Upon activating the geolocation tool on a raster layer, a small dialog will appear floating on the view with a series of text entries and a button bar. Also the raster layer that has been selected acquires the capacity to be moved, rotated or moved by clicking and dragging it to the right place on the view.

Transformation on the view

geolocalizacion.img/es/geolocalizacionRaster.png

Vista Panel Image Geolocation

Transformation by introducing coordinates

From the geolocation floating box we can modify the transformation of the image on the view. The text entries marked as X, Y, Pix X, Pix Y, Rot X, Rot Y contain the current position in upper left X coordinate, upper left Y coordinate, pixel size in X, pixel size in Y, rotation in X and rotation in Y respectively. If we change these values from the keyboard these will be updated on the image by pressing the "Enter" key while the cursor is inside the text box or when it loses focus.

Controls of the geolocation box

Centering the raster layer in the view: With this button we center the raster in the current view regardless of zoom it has selected. The coordinates will be automatically calculated.


Georeferencing
Description

Direct traslation from Google traslator

To launch the georeferencing dialog it is used the dropdown toolbar selecting the "Geographic Transformations" button on the left and "Georeferencing" from the dropdown button on the right.

georreferenciacion.img/es/barraGeorref.png

Georeferencing Tool

Initially we must decide what type of georeferencing to implement, "reference maps" or "without reference maps".

Georeferencing with "Mapping Reference"

Start Dialogue georeferencing

To implement this type of georeferencing is imperative that we have previously charged in a view mapping that we will provide a geographic reference for taking control points. In case of not having it will close the options dialog georeferencing and proceed to prepare for the hearing. Once we have the view with reference maps georeferencing tool launched will see that the option "reference maps" is checked by default. Below is a dropdown menu which lists the views that gvSIG has at that time. If you have several it must select a view which is our base mapping for decision-points.

georreferenciacion.img/es/Launcher.png

Dialog Georeferencing

In the panel marked "a georeferenced file" pops up a dialog for selecting the file for which you want to create checkpoints and later georeferencing.

The panel labeled "Output File" we must put the path and file name destination if the georeferencing is done with resampling. This option can vary from box options once we are inside the application, so it is not essential to a correct value at the moment, but if must be done before the end of the process.

The panel "georeferencing algorithm" select how we will get the output result. There are two possibilities, "affine transformation" and "polynomial transformation".

The affine transformation applied to raster an affine transformation only to the calculations performed with the control points taken. The affine transformation applied will be allocated "on the fly" for the display and the output image is the same as the input. The result of this transformation is therefore a georeferencing file. Keep in mind that this type of transformation is limited and the user will be responsible for selecting the most convenient transformation in each case.

The polynomial transformation involves a resampling of the input image taking into account the reference control points and obtaining an output image with deformations necessary to adapt to the new location. If you select this option we will be forced to decide the degree of transformation that we apply and the type of interpolation that we want to apply for calculating new pixels. Depending on whether you choose one degree or another need a minimum number of control points for them. This number of points required is given by the formula (order + 1) * (order + 2) / 2, ie for a polynomial of degree one will be needed at least three points, to grade two will need six points for third grade ten points ... The interpolation method affects the way we calculate the information that we have not. When an image georeferenced output image has deformations with respect to the original there are areas where no information is available. These can not be empty with what must be calculated from the areas where we know. These calculations can be performed by various methods, the simplest of these is "Nearest neighbor" which will be unknown pixel information closest known pixel. Other methods such as "bilinear" or "bicubic" make calculations using the known group of pixels surrounding the unknown. These other methods give a more relaxed but it is slower in its implementation. This option can vary from box options once we're inside the application.

The panel "Pixel Pitch" is the pixel size information of the output image. In principle this will be calculated from the input image but can be changed manually. This option can vary from box options once we are inside the application, so it is not essential to a correct value at this time.

The views

Executing the application are two views. The left contains the base mapping that we carry in the gvSIG view of the right and the image we want to georeference. Both have a control bar on the right for view actions. Also in the upper left corner are the coordinates of the mouse cursor. In reference mapping coordinates are those of the real world. In the image to be georeferenced coordinates in pixel coordinate on the upper left.

Cursor Zoom

In the central part appears a cursor with a central window. The window cursor is active when the view can be resized and moved. The contents of this window will be on display in the zoom windows. Ca da vista has its associated zoom window at the bottom. Par resize window cursor select the view you want by clicking on it then bring the mouse to the edges of the window until the pointer changes to horizontal or vertical arrows. Now we click and drag to force the resizing. To move the cursor window select the view you want by clicking on it then bring the mouse to the corners of the window until the pointer changes by crossed arrows. We must now drag and drop to force displacement.

View Controls

There are six controls to handle the zoom level and position of the view mapping Increase the level of zoom: the zoom level increases by multiplying by 2 the current level.

georreferenciacion.img/es/mas.png

Decrease zoom level: it decreases the zoom level by dividing by 2 the current level.

georreferenciacion.img/es/menos.png

Zoom area selection: Activates a tool on the hearing in order to make a rectangle the area we want to see enlarged.

georreferenciacion.img/es/zoomSeleccionArea.png

Full Zoom: Put a zoom level so that you can view the entire mapping.

georreferenciacion.img/es/zoomTodo.png

Zoom Previous: Sets the zoom level that you previously selected.

georreferenciacion.img/es/zoomPrev.png

Displacement: clicking and dragging on the scroll view mapping.

georreferenciacion.img/es/pan.png

Zoom controls

Each view has an associated georeferencing zoom window centered over the cursor. When we move the cursor on the sale of view varies the position where the zoom and focus when we change the window size changes the zoom level. In the upper left corner of the window coordinates of the mouse cursor as in the overview.

georreferenciacion.img/es/zoom.png

Zoom box associated with the views

Checkpoints

A control point is an entity that provides a correspondence between a geographic coordinate and pixel coordinate. Control points are represented in raster geographic view as Blue-and red circles respectively. To add a new control point is selecting "New" in the table control. This makes a new entry in the table appears. A control point is associated with a table entry. By selecting "New" automatically creates a point at coordinates 0, 0 for both views and will activate the tool "move point". Now clicking on the view point where we will move puncture. We assign the coordinate point numerically by writing directly on the input value in the table (X for the geographic coordinates X, Y geographic coordinate for Y, X 'for X and Y pixel coordinate' for the pixel Y coordinate). The points can also be moved by clicking and dragging on them. This may be done both in hearings and in the zooms.

georreferenciacion.img/es/Georreferenciacion.png

The process of georeferencing. Sights and points of control

Options

The panel labeled "Output File" we must put the path and file name destination if the georeferencing is done with resampling.

The panel georeferencing algorithm "select how we will get the output result. There are two possibilities, "affine transformation" and "polynomial transformation".

The affine transformation applied to raster an affine transformation only to the calculations performed with the control points taken. The affine transformation applied will be allocated on the fly for the display and the output image is the same as the input. The result of this transformation is therefore a georeferencing file. Keep in mind that this type of transformation is limited and the user will be responsible for selecting the most convenient transformation in each case.

The polynomial transformation involves a resampling of the input image taking into account the reference control points and obtaining an output image with deformations necessary to adapt to the new location. If you select this option we will be forced to decide the degree of transformation that we apply and the type of interpolation that we want to apply for calculating new pixels. Depending on whether you choose one degree or another need a minimum number of control points for them. This number of points required is given by the formula (order + 1) * (order + 2) / 2, ie for a polynomial of degree one will be needed at least three points, to grade two will need six points for third grade ten points ... The interpolation method affects the way we calculate the information that we have not. When an image georeferenced output image has deformations with respect to the original there are areas where no information is available. These can not be empty with what must be calculated from the areas where we know. These calculations can be performed by various methods, the simplest of these is "Nearest neighbor" which will be unknown pixel information closest known pixel. Other methods such as "bilinear" or "bicubic" make calculations using the known group of pixels surrounding the unknown. These other methods give a more relaxed but it is slower in its implementation.

The panel "Pixel Pitch" is the pixel size information of the output image. In principle this will be calculated from the input image but can be changed manually.

The panel labeled "Options" contains settings of a different nature. Since we can change the background color of view, the text color of the views. The "show the number of graphically checkpoint" will be displayed or hidden by the control point a point that indicates the corresponding point number. "Add the CSV file errors" will be generated when this type of text files with all the control points we can ignore the file or add the calculated errors. The "Focus the selected point view" makes automatically every time we select a point on the table the view is focused on this. The effect is much as if the tool center point was always active. The "error threshold for the warning," assigns the value at which the error appears in red on the table.

georreferenciacion.img/es/Opciones.png

Options for georeferencing

Points Table

The points table is below the sights and initially will be empty. Each table entry corresponds to a checkpoint. It appears all the information related to a point. This table can see it folded its default state or maximized. In its maximized state are folded more information. On the left side of the row there is a check to activate and deactivate the current row. This means that this point will not be displayed graphically or be taken into account for calculation errors and will be prosecuted to do a test. The information can be found in the points table on each point:

  • Number of point
  • Real coordinate X
  • Real coordinate y
  • Coordinate pixel X
  • Coordinate pixel Y
  • Error in X
  • Error in Y
  • Total RMS error for that point

The quality of the geometric correction can be estimated based on the mean square error RMS error and the contribution of each point. When the contribution to RMS of a point is high, this may indicate that the correspondence of points was poorly selected and the point is not well suited to model transformation between image and map or other information used as reference. The points with high contribution that exceeds a certain threshold can be deleted or deactivated, and calculate the RMS. However, when we are fully confident of the location of a point, and to find you, the RMS is triggered, it may be possible that the geometric model does not resolve the local arrangements, for which they may need a better model, which means, put some more points, right on the problem area.

There is also a global RMS error in an external text field for all points.

georreferenciacion.img/es/TablaPuntos.png

Control Panel points

Controls

Tool center point: When you press the focus control to the view point that is selected.

georreferenciacion.img/es/CentrarPunto.png

Georeferencing operation completes. Before you ask if we carry on the gvSIG view the results of the last trial. You'll also want confirmation of application output.

georreferenciacion.img/es/FinalizarGeoref.png

Launching the options dialog.

georreferenciacion.img/es/OpcionesBoton.png

Make a test with the control points currently entered. If there are not enough for the specified algorithm will warn. The result is that applying the transformation and loading the transformed image on the view with the reference maps.

georreferenciacion.img/es/Procesar.png

Save the control points in the metadata file attachment with the raster.

georreferenciacion.img/es/SalvarRMF.png

Retrieves the control points that are in the metadata file attached to raster.

georreferenciacion.img/es/CargarRMF.png

Ends the test of processing the raster. Eliminate the test image loaded in the view with the mapping.

georreferenciacion.img/es/TerminarTest.png

When the button "Select point" we are active, clicking on the view assigning the selected point on the table at that time to the position.

georreferenciacion.img/es/SeleccionarPunto.png

Sequence capture control points

There may be ways to capture control points with the tools available. An example would be the following sequence of actions:

  • Click "New" in the table of control points. This will create a new row is selected in the table. In addition the tool "Move Point" is selected.
  • Click with your mouse pointer over the view to locate the point raster.
  • Click with your mouse pointer over the view with reference maps to locate the point.
  • Push the button "Refocused selected view point" to place the checkpoint in the center and appears in the zoom window.
  • With the tool of choice for area Zoom "or" increase the level of zoom "or" Decrease the zoom level "we can set the desired zoom level until the controls of" Zoom "we have an optimal resolution level approximate.
  • Click and drag the control point in the zoom window to place it more precisely. The accuracy depend on how correct is selected previous zoom level.
  • Use the zoom tools to return the view to a wider zoom level and to allow a new control point.
  • To return to one point and reset the Selection Click on the row of the table, click "Center views the selected point, adjust the zoom level to zoom tools and we'll move the view by clicking and dragging on the window zoom for greater accuracy.

Georeferencing with resampling

Two types of processing for raster. If selected in the options the affine transformation the image obtained is not wide and applies an affine transformation on the view. This transformation is a scaling, displacement, rotation and deformation in the direction of axis X and / or Y axis The transformation with resampling involves generating a new image from the original on which areas can appear empty. These areas are due to the fact that the resulting image should be rectangular but the area covered by the data processing may not have applied this same way.

georreferenciacion.img/es/GeorreferenciacionRemuestreo.png

Results georeferenced image with resampling

Once the process of georeferencing the raster generated and loaded in the view we can apply a transparency per pixel to eliminate the empty areas.

georreferenciacion.img/es/GeorreferenciacionRemuestreoConTransp.png

Image georeferenced, with application of transparency

Georeferencing without "Mapping Reference"

The georeferencing without reference maps is useful when you do not have imagery that guide us to assign the control points. We will have to allocate the actual coordinates directly by typing its value. In this case it is useful in view of the left so it will allow more space for the raster and the points table. The operation is very similar to the two views just that when you select the point on the reference maps have to type the entry of the table directly.

The operation of other controls is the same as with reference maps.


Modifying georeferencing upon loading

When you load a file that is not georeferenced, gvSIG can prompt you to enter the coordinates manually. For this the option has to be activated in preferences, which is disabled by default. The option that needs to be activated is "Pedir las coordenadas al cargar un raster que no tiene georreferenciación".

georreferenciacion.img/es/PreferenciasModificargeorreferenciacioncargacapa

Raster Preferences. Pedir las coordenadas al cargar un raster que no tiene georreferenciación

In this case we will see a dialog with the message: "Name of layer. The layer has no georeferencing. Do you want to enter them manually?". If you answered "No", the load is carried out with the coordinates (0, 0) in the top left corner (width in pixels, height in pixels) in the lower right corner. If the answer is "Yes", then a dialog to enter the coordinates of the raster will show. In this one must be careful to enter valid data to avoid erroneous results. The dialog has two tabs from which we can enter the coordinates in the form of affine transformation or the upper-left and bottom-right corners. In the first mode will need the X and Y coordinates of the upper right corner of the original raster, the pixel size in X, the pixel size on Y, X rotation and the rotation Y.

georreferenciacion.img/es/modifTransformacionRaster.png

Loading raster without georeference. Transformation

In the second mode we only introduce the coordinates of the corners in the order indicated by the graphs.

georreferenciacion.img/es/modifEsquinasRaster.png

Loading raster without georeference. Corners


Image reprojection

For the reprojection of raster layers, gvSIG uses the GDAL library. The reprojection process can be launched in two different ways: By activating the reprojection icon from the raster toolbar for images that have already been loaded to the view, or by reprojecting the layer before it is loaded to the view if this is needed.

The GDAL library does not support ecw, mrsid or jpeg2000 images and therefore images in these formats cannot be reprojected.

To launch the reprojection dialog from the raster toolbar, select "Geographic transformations" on the left drop-down button and "Reproject layer" on the drop-down button on the right. Make sure that the layer that you want to reproject is set as the current layer in the text box.

reproyeccion.img/es/barra_reproyect.png

Image reprojection icon on the raster toolbar

When launching the reprojection function from the raster toolbar, a dialog opens which shows the projection information of the input image as "source projection". The source projection cannot be changed as it is assumed that the input layer has been loaded into the view with the correct projection. Under "target projection" the projection of the output image can be set by the user through the standard gvSIG dialog for CRS and transformations. It should be noted that not all transformations are supported; the projection options depend on the GDAL reprojection library.

The output layer can be saved on disk or opened in memory as temporary file. When the first option (which is the default) is selected, the user is prompted for a file name and path. Then, the reprojection process starts and when this is finished, it will ask whether you want to add the new layer to the TOC.

reproyeccion.img/es/reproy.png

Reprojection dialog

NB: When reprojecting an image, the used transformation is "EPSG Transformation"; with raster layers the other transformations (manual, composed or grid) can not be used.

Images can also be reprojected before loading them to the view. To do this, you will need to have the option "Ask for projection when the raster loaded has different projection from view", located in the raster options section of the Preferences dialog, selected. If this option is selected and a raster with a different projection than the view is loaded, a dialog is opened with projection options. The default option is to load the layer while ignoring the projection, but you can reproject the layer by selecting the option "Reproject raster to the view projection". Then, the same Reprojection dialog is shown, but in this case the "target projection" is fixed to the projection of the view, and the "source projection" can be changed, as in some cases the projection of an image may not have been set correctly or the needed projection information maybe missing.

After accepting the settings, the reprojection process will start and the layer will be added to the TOC.

reproyeccion.img/es/reproy_carga_capa.png

Reprojection options when loading an image to the view


Cached time 11/22/13 02:05:33 Clear cache and reload